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Predicted strain-induced isosymmetric ferri-to-ferroelectric transition with large 

piezoelectricity 
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We identify a first-order, isosymmetric, transition between a ferrielectric (FiE) and ferroelectric 
(FE) state in A-site ordered LaScOs/BiSc03 and LalnOs/BilnOs superlattices. This unusual 
electrical analogue to the pervasive ferrimagnetic state in magnetic oxides relies on a geometric field 
imposed by octahedral rotations and enabled by cation ordering — it decouples the local electric 
dipoles on the cationic sublattices and directs their cooperative alignment. Epitaxial strains less than 
2% are predicted to be sufficient to transverse the phase boundary, across which we capture a ~ 5X 
increase in electric polarization. We describe the electromechanical responses and show that at fixed 
composition the FiE-to-FE transition may be superior to the polarization reversal found in Pb-based 
perovskite oxide solid-solutions, enabling piezoelectric technologies under extreme conditions. 

PACS numbers: 77.84.-s, 81.05.Zx, 71.15.Mb, 68.55.-a 



Introduction. — Ferroic materials are characterized by 
the loss of spacial and/or temporal symmetries. Inde- 
pendent of the microscopic ordering, symmetry removal 
provides for functional responses that couple multiple 
degrees of freedom, activating new physical tensor proper- 
ties. Perovskite ABO% oxides with broken space inversion 
symmetry often exhibit a spontaneous and switchable 
electric polarization: Ferroelectric (FE) ordering induces 
an alignment of local electric dipoles, in analogy to the 
cooperative ordering of local magnetic moments in ferro- 
magnetic crystals without time-reversal symmetry. Much 
rarer in crystalline materials with an electric polarization, 
however, is the appearance of a ferrielectric (FiE) state, 
vis-d-vis ferrimagnetism, where local electric dipoles of 
different magnitude are anti-aligned to yield a net non- 
zero electric polarization [1-4]. The underlying reason 
being that the long-range Coulomb forces in oxide-based 
dielectrics favor the cooperative alignment of all electric 
dipoles in the crystal [5] through cation displacements 
that occur against an oxygen ligand framework [6]. 

The recent discovery of a class of hybrid improper 
ferroelectrics (HIF) [7, 8], whereby the electric polariza- 
tion originates from a non-zero cancellation of electric 
layer polarizations, or unequal local electric dipoles mo- 
ments, p, among three cation sublattices is a realization of 
the FiE state [9]. In cation ordered (A^A f )B 2 OQ oxides, 
equivalent ly 1/1 superlattices constructed interleaving 
ABOs/A'BOs blocks along [001], the layer dipoles are 
restricted by symmetry to be directed along the pseudo- 
cubic (110) directions. At present, all known examples 
of this type are found to have two of the three possi- 
ble cation-sublattice dipoles aligned with each other, e.g. 
[110], while the third is anti-parallel to the others, [110] 
[10]. A rotocrystalline anisotropy stabilizes the FiE state, 
because the electric dipoles of the A-site sublattice are 
constrained by the "sense" of the orthorhombic a~a~c + 
BOq rotations [11] in Glazer notation [12]. The electric 



polarizations of each sublattice, therefore, may be tuned 
independently — the A and B cation sublattice polariza- 
tions rely on the geometric field imposed by the rotations, 
which can decouple the electrostatic alignment of the 
layer dipoles. Can a suitable "field" convert the FiE state 
directly into a FE state without passing through a para- 
electric intermediate? If answered in the affirmative, then 
a rotation-enabled route to achieve a significant polar- 
ization enhancement and large piezoelectric response in 
synthetically ordered perovskites would be realized. 

In this Letter, we use first-principles density functional 
calculations to demonstrate that the polarization switch- 
ing between a FiE and FE phases in the (A 1 A')B20q 
perovskites is accessible with substrate-induced epitaxial 
strain. Specifically, we uncover a first-order isosymmet- 
ric phase transition and evaluate the feasibility of the 
FiE-to-FE transition through a minimum energy path 
(MEP) trajectory. We find enhanced electromechanical 
responses when passing through the metastable intermedi- 
ate structure, which resemble and are comparable to the 
the composition-dependent morphotropic phase bound- 
aries (MPB) in solid-solution PbZri_ x Ti^0 3 [13, 14], 
where a large piezoelectric response is driven by rota- 
tion of polarization along an easy rotation path between 
polar rhombohedral and tetragonal phases through a mon- 
oclinic intermediate [15, 16]. Lastly, we describe how the 
transition manifests as an effective MPB without requir- 
ing a change in octahedral rotation sense. It has the 
practical advantage of occurring within a single phase 
material without a change in symmetry and at integral 
stoichiometry — features useful for facile switching and 
mitigating piezoelectric fatigue [17]. 

Models and Simulation Details. — We select ultra-short 
LaMOa/BiMOs perovskite superlattices, with trivalent 
cations M = Ga, Sc and In, respectively, to guarantee 
the insulating behavior required to sustain an electric 
polarization. The 6s 2 stereochemical lone pair of Bi pro- 
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FIG. 1. (Color online) (a) The Pmma paraelectric structure 
of the LaScOa/BiScOs superlattice with the a~a~c° rotation 
pattern (a). The displacement patterns of the corresponding 
soft phonon modes (b): cji with anti-parallel La-O/Bi-0 in- 
plane displacements (left), and uju with parallel La-O/Bi-0 
in-plane displacements (right). 



motes structural distortions and enhanced polarizability 
[18]. Lastly, each bulk A 3+ M 3+ 3 oxide forming the 
superlattices exhibits the targeted orthorhombic a~a~c + 
octahedral rotation pattern (Fig. SI of Ref . 19) conducive 
toHIF[ll]. 

Our density functional theory (DFT) calculations are 
performed with the QUANTUM ESPRESSO code [20], with 
optimized nonlocal norm-conserving [21] pseudopotentials, 
a 60 Ry plane-wave cutoff, and the PBEsol functional 
[23]. The Brillouin zone (BZ) integrations use a 8x8x8 
Monkhorst-Pack /c-point grid [22]. Low-energy phases are 
obtained by performing complete structural relaxations 
on soft-mode generated structures until the Hellmann- 
Feynman forces on the atoms are less than 0.1 meV A -1 . 
For the strained films, the out-of-plane lattice constant 
is relaxed until the stresses are less than 0.1 kbar. The 
electronic contribution to the polarization is calculated 
following the Berry phase formalism [26]. The minimum 
energy path (MEP) between the stable phases is simulated 
using the nudged elastic band (NEB) method [27]. 

Ground State Structures. — We find the lowest energy 
structures by computing the phonons in the BZ of the 
parent PA/mmm phase, i.e., the distortion- free structure 
obtained from A/ A' cation ordering along [001]. We find 
four unstable vibrational phonons in each superlattice: 
They are characterized as either in-phase (irreducible 
representation M% relative to Pm3m) or out-of-phase 
rotations (R^) of M0 6 octahedra, or polar inversion 
symmetry lifting cation displacements along [110] (TJ ) 
and [001] against the oxygen framework, respectively [19]. 
For all superlattices, the most unstable modes are out-of- 
phase rotations, a~a~c°, with the mode- frequency of the 
gallate smaller than the scandate or indiumate due to the 
larger tolerance factor of the former [9] . 

We next freeze-in each soft mode into the P4/mmm 
superlattices followed by full structural relaxations to 
obtain lower symmetry distorted phases. For each su- 
perlattice, we find a stable non-polar Pmma phase with 
the a~a~c° rotation [Fig. 1(a)]. We now take this as the 



TABLE I. The relative energetic stability (in meV/f.u.) with 
respect to the ground state Pmc2i (phase I) structure for 
various low-energy phases specified by Glazer tilt; subscripts 
indicate the direction of the electric polarization following the 
convention in Ref 28. The average tolerance factor (r) for 
each perovskite superlattice is given in parentheses. 





Pmma 


[Pmc2i]i 


[Pmc2i] n 


Pc 


System / Rotation 


a a Q c 


a7_a7_c+ 


a+a~c+ 


a7a7c7 


Laln0 3 /Biln0 3 (0.87) 
LaSc0 3 /BiSc0 3 (0.89) 
LaGa0 3 /BiGa0 3 (0.94) 


437 
338 
97 







109 

57 


179 
170 



relevant paraelectric phase for which transitions between 
lower-symmetry polar phases could occur. We find two 
unstable vibrational modes in the Pmma superlattices 
(Table SII in Ref. 19): uu\ and uju correspond to atomic 
displacement vectors that are either a combination of 
in-phase octahedral rotations and anti-parallel A cation 
displacements, or cooperative polar cation displacements 
moving in the direction opposite of the oxygen anions 
[Fig. 1(b)], respectively. 

We now combine the Pmma structures with the dis- 
placement vectors for uo\ and uj\\. The symmetry of the 
superlattices all reduce to the same polar space group 
Pmc2\ (no. 26), where the polar axis is along the pseu- 
docubic [110] direction, despite the physical difference 
in the displacement vector. Hereafter, we distinguish 
between the two symmetry equivalent Pmc2\ phases as 
(I) or (II), appended as a subscript to the space group 
label. After performing a full structural relaxation, we 
find that the [Pmc2\]i phase is the global ground state for 
all superlattices (Table I). Except for LaGaOs/BiGaOa, 
[Pmc2i]u is the next lowest energy phase, e.g., 57 meV 
per formula unit (f.u.) higher in energy than [Pmc2i]\ 
in LaScOa/BiScOs [29]. The small energy barrier sug- 
gests transition between the two Pmc2\ phases may be 
achieved by thin film epitaxy. 

To understand the microscopic origin of the energy dif- 
ference between the two isosymmetric Pmc2\ phases, we 
focus on the LaScOs/BiScOa superlattice, and examine 
its local atomic structures in more detail. Fig. 2(a) and 
(b) shows both polar phases projected along the [100] and 
[011] directions. Each possesses the same a~a~c + rota- 
tions, but different cation displacement patterns, which 
become readily apparent in the layer-resolved electric 
polarizations [Fig. 2(c)]. [Pmc2i]i is described by a ferri- 
electric state, whereby a non-cancellation of anti-aligned 
layer electric dipoles from the alternating LaO or BiO 
planes produces a overall non-zero electric polarization. 
In [Pmc2i]n, PBiO switches from anti-parallel to parallel 
alignment with PLaO- All cations displace in the same 
direction, as a conventional ferroelectric, to give a >5X 
polarization enhancement. 

Phase I is obtained by condensation of the ljj soft mode, 
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FIG. 2. (Color online) Atomic structures of the two isosym- 
metric polar ground states (cl+cl+Cq) with Pmc2\ symme- 
try: phase I (upper) and phase II (lower). Both top view 
(a) and side view (b) are shown to emphasize cation dis- 
placement directions (arrows), (c) Layer-revolved polarization 
Player = ^ _1 J2i u i^i^ where Q is the cell volume, m is the 
i cation displacement and Z* is the Born effective charge. 
Phase I is ferrielectric (total polarization, P — 0.12 C/m 2 ) 
with anti-parallel pLaO and pbio, while Phase II is a normal 
ferroelectric (P — 0.63 C/m 2 ), where both pLaO and pBio are 
oriented in the same direction. 



and the polar displacements are stabilized from a combina- 
tion of the a~a~c° and a°a°c + octahedral rotations (Fig. 
SI). Interestingly, the extraordinary large energy gain 
between polar and paraelectric states of LaScOs/BiScOa 
(338 vs. 48 meV/f.u. as in BaTi0 3 [30]) indicates the 
polarization should persist to high temperatures, similar 
to what found in other HIF oxides [31]. The polarization 
in phase II arises from the polar instability ujjj in Pmma, 
where cooperative displacement of all cations occur in 
the same direction and are opposite to that of the oxygen 
atoms, as found in proper ferroelectrics, e.g., BaTiOs. 
Moreover, phase II has a suppressed in-phase rotation 
amplitude compared to phase I (Fig. S2) in the presence 
of a greater electric polarization, indicating the trilinear 
coupling term connecting the octahedral rotations to the 
polarization contributes less to the ground state [8, 9]. 

Strain-induced FiE-to-FE Transition. — We calculate 
the evolution of the total energy for the Pmc2\ phases 
with respect to applied epitaxial strain [Fig. 3(a)]. We 
find in LaScOs/BiScOa, the transition from [Pmc2\]i to 
[Pmc2i]u occurs approximately at 1.0% tensile strain, 
accompanied by an abrupt change in the overall polar- 
ization [Fig. 3(b)]. The transition between the Pmc2\ 
phases is first-order and isosymmetric [32]. Critically, 
no octahedral rotation-reversal is necessary; only a re- 
duction in amplitude is required, indicating polarization 
reversal between the two states should be experimentally 
simpler to achieve than in other HIF. [8] (See Fig. S2 for 
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FIG. 3. (Color online) Strain-dependent phase stability and 
electronic properties of LaScOs/BiScOs: (a) total energies for 
the paraelectric and polar structures, (b) total polarization 
[the bold black line transverses the polarization of the lowest 
energy structures in (a)], (c) static dieletric tensor dj (e xx —e yy , 
filled circles), and (d) the piezoelectric stress constants eij 
(e33, filled circles; e^i empty circles) for Pmc2\ phase I (red) 
and II (blue). The isosymmetric transition occurs at 1.0% 
tensile strain (broken lines). 



detailed structural analysis.) The ferri-to- ferroelectric 
transition also occurs in the indium superlattice, albeit 
at a larger strain (~2.0%), yet it is not found in the 
gallium compound [19]. The larger tolerance factor for 
the LaGaOa/BiGaOa reduces the energy gain between 
Pmma and [Pmc2i]i (Table I), and makes the [Pmc2i]u 
phase unstable. In the remainder of this Letter, we fo- 
cus on the strained LaScOs/BiScOs with an emphasis 
on response property enhancement near the FiE-to-FE 
transition. 

The evolution of the response functions, including the 
dielectric permittivity and piezoelectric response, as a 
function of strain are computed at the DFT-level [33]. 
Fig. 3(c) shows that the in-plane dielectric susceptibility 
e xx and e yy are maximum at the FiE-to-FE transition 
at 1.0 % tensile strain in both Pmc2\ phases, consistent 
with a first-order transition at this critical strain [34]. In 
Fig. 3(d), we also plot piezoelectric stress constants 633 
and ei3, which evolve smoothly with respect to strain. 
The non-zero piezoelectric strain coefficients dij for the 
FiE and FE phases at 1.0% strain are also reported in 
Table II. Compared to proper FE perovskites, each polar 
LaScOs/BiScOs displays a modest piezoelectric response. 
In polar Pmc2\ state where the polarization is induced 



TABLE II. Calculated non-zero piezoelectric stress constants 
Cij (C/m 2 ) and piezoelectric strain coefficients dij (pC/N) for 



each phase 


of LaScO: 


3/BiSc0 3 


at 1.0% tensile strain. 








[Pmc2i]i 


[Pmc2 1 ] u 




index (ij) 


Cij 




dij eij 


dij 



31 
32 
33 

24 
15 



0.50 
1.05 
1.34 
0.93 
1.10 



-1.05 
2.33 
5.23 
14.57 
14.30 



-0.19 
0.77 
2.45 
1.03 
1.23 



-6.16 
0.91 
15.66 
17.96 
21.61 



by octahedral rotations, the dominant strain-rotation cou- 
pling will mitigate the effective coupling between strain 
and polar displacements [35, 36], making hybrid-improper 
ferroelectric oxides weak piezoelectrics for transducer ap- 
plications [31]. Despite the small piezoelectric moduli 
found in the individual [Pmc2\]i and [Pmc2i]u phases, 
one should also consider enhancements due to the FiE-to- 
FE transition in describing the macroscopic functionality. 

Piezoelectricity at the FiE-to-FE Transition. — To eval- 
uate the consequence of the FiE-to-FE transition on 
piezoelectricity, we first map out the energy landscape 
in the vicinity of the low energy phases (Fig. 4). The 
two-dimensional energy contours at 1.0% strain for 
LaScOa/BiScOs are computed by varying the amplitude 
of the eigenvectors of the unstable Pmma phonons. We 
find that the FiE [Pmc2\]i and FE [Pmc2i]u phases are 
located at minima near the major axes and are separated 
by a energy barrier. 

We estimate the feasibility of the FiE-to-FE transition 
by calculating the energy barrier separating the two stable 
structures along a MEP obtained from the NEB method 
[Fig. 4(b)]. We find a small, nearly strain-independent 106 
meV/f.u. barrier separating the phases [37], and compa- 
rable to the polarization rotation barrier in piezoelectric 
PbTiOa [38], suggesting the FiE-to-FE transition should 
be energetically accessible [39]. 

We emphasize that the polarization rotation at the 
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FIG. 4. (Color online) The calculated two-dimensional energy 
surface contours for LaScOs/BiScOs with respect to the soft 
phonons of Pmma (a); energy contours in meV. The MEP 
between two polar states for LaScOs/BiScOs under tensile 
strain (b) reveals a small energy barrier. 
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FIG. 5. (Color online) The calculated (a) total polarization, 
(b) layer-revolved polarizations and (c) major crystallographic 
mode amplitudes for each image along the MEP connecting 
phase I and II in LaScOs/BiScOs at 1.0 % strain. The position 
of the saddle point along the MEP is shown as dashed line. 
Electronic-localization-function (ELF) surfaces within the BiO 
layer for phase I, the intermediate structure at the saddle 
point, and phase II (d). Across the FiE-to-Fe transition, the 
Bi electron lone pair (small lobe, arrowed) alignment reverses. 



[Pmc2i]i — >• [Pmc2i]u transition [Fig. 5(a)] along the 
MEP is due to a change in direction of Bi displacements 
along [110] direction in (001) planes above and below 
the Sc02 layers (Fig. 2). peio switches from anti-parallel 
to parallel alignment with pLaO along the polar axis by 
decreasing its amplitude in the (001) plane until reach- 
ing the saddle point in the energy landscape [Fig. 5(b)]. 
Here, a nearly random distribution of the Bi lone pairs 
within the planes emerges [Fig. 5(d)], and we find a trans- 
formation of BiOg polyhedra into highly distorted BiC^ 
tetrahedra with enhanced Bi 3+ electron lone pair activity 
in [Pmc2i]u [Fig. S4 and Fig. S5]. After the saddle point, 
PBiO continues to increase while pLaO remains nearly un- 
changed. 

The change in mode amplitudes along the MEP is dom- 
inated by a large enhancement of T^ polar displacement 
(>5X), with the moderate change in octahedral rotations 
[~25%, Fig. 5(c)]. This structure evolution indicates 
the lattice strain-polarization coupling can exceed the 
strain-rotation interactions across the FiE-to-FE transi- 
tion, despite reduced coupling in each individual phase. 
Because the polarization remains finite - the FiE state is 
directly converted into FE state without passing through 
any paraelectric intermediate - the changes in polyhedral 
coordination at the transition can produce a large induced 
electric polarization upon application of a stress. The 



FiE— ^FE transition occurs through an easy polarization 
switching path with a large polarization enhancement and 
strong strain-polarization coupling. Indeed, the piezoelec- 
tric moduli are enhanced from 3-6X across the transition 
(Table II), which is qualitatively similar to conventional 
MPB oxides [14], yet occurs at fixed composition. 

Outlook & Summary. — We used first-principles calcula- 
tions combined with group theoretical studies to uncover 
a first-order, isosymmetric, structural transition between 
a FiE and FE state in artificial A-site ordered perovskite 
superlattices. For LaScOa/BiScOs, we predict 1% tensile 
epitaxial strain is sufficient to induce the FiE-to-FE tran- 
sition with electrical features that could enable rational 
technology design: A single phase material with a two- 
state polarization enabled by strong strain-polarization 
coupling, and simultaneously large piezoelectric responses 
driven by reversal of cation displacements without changes 
to the sense of octahedral rotations or crystallographic 
symmetry — features absent in conventional ferroelectric 
ceramics with MPB. Second, owing to the extraordinary 
large energy gain between the polar and paraelectric 
states, the FE-to-FiE transition should persist up to high 
temperature (T > 500° C), where Pb-based piezoelectrics 
responses are impaired. Finally, ferrielectric materials 
should exhibit similar electrical switching behavior as an- 
tiferroelectrics, e.g., double hysteresis loops that improve 
energy density capacities [40], yet with the advantage 
of polarization reversal at much lower fields, enabling 
practical integration for energy harvesting devices. 
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